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Abstract 

I We present a theoretical prediction for the photon spectrum in radiative T decay. Parts of 

' the spectrum have aheady been understood, but an understanding of the endpoint region has 

^ remained elusive. In this paper we provide the missing piece, and resolve a controversy in the 

Q . literature. We treat the endpoint region of T ^ Xj decay within the framework of the soft- 

Q , collinear effective theory (SCET). Within this approach the T structure function arises naturally, 

l/^ I and kinematic logarithms are summed by running operators using renormalization group equations. 

In a previous paper we studied the color-octet contribution to the decay. Here we treat the color- 

^ ■ singlet contribution. We combine our result with previous results to obtain the T X'j spectrum. 

■ We find that resumming the color-singlet contribution in the endpoint gives a result that is in much 

, better agreement with the data than the leading order prediction. 
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I. INTRODUCTION 



The large mass of the constituent quarks in quarkonium makes this system simple enough 
to use as a testing ground for theory. For example, the spectrum of quarkonium has been 
extensively studied in Lattice QCD, allowing qualitative investigations of systematic errors 
and extractions of parameters (such as the strong coupling constant and the b quark mass). 
In addition, the gross under-estimate of the production of tp' at the Tevatron helped establish 
the validity of a Non-Relativistic effective field theory of QCD (NRQCD) [0, 0. While much 
has been learned about the charmonium and bottomonium systems, it is still a useful probe 
for theory. 

Inclusive decays of quarkonium are understood in the framework of the operator product 
expansion (OPE), with power-counting rules given by NRQCD. The OPE for the direct 
photon spectrum of T decay is [Q 

^ = ^C„(M,^)(T|a|T), (1) 

n 

where z = 2E^/M, with M = 2mh. The Cj are short- distance Wilson coefficients which 
can be calculated as a perturbative series in as{M), and the O are NRQCD operators. 
NRQCD power counting assigns a power of the relative velocity, f <^ 1, of the heavy quarks 
to each operator, which organizes the sum into a power series in v. At leading order in v 
only one term in the sum must be kept, the so called color-singlet contribution. The color- 
singlet operator Oi{^Si) creates and annihilates a quark-antiquark pair in a color-singlet 
^Si configuration, and is multiplied by the color-singlet Wilson coefficient, which at leading 
order is proportional to a^(M). 

This simple picture of the photon spectrum in inclusive T decays is only valid in the 
intermediate range of the photon energy spectrum (0.3 ^ z ^ 0.7). In the lower range, 
z ^ 0.3, photon-fragmentation contributions are important ||^, ^. At large values of the 
photon energy, z ^ 0.7, both the perturbative expansion and the OPE ^ break down. 

The breakdown of the OPE and the perturbative expansion is a consequence of NRQCD 
not containing the correct low energy degrees of freedom to describe the endpoint of the 
photon spectrum. The effective theory which correctly describes this kinematic regime is 
a combination of NRQCD for the heavy degrees of freedom, and the soft-collinear effective 
theory (SCET) P, |^, H, |^ for the light degrees of freedom. In a previous paper |jTO| we applied 



SCET to the color-octet contributions to radiative T decay. First we used the SCET power 
counting to show that in the endpoint region there is a color-octet contribution which is 
the same order as the color-singlet contribution. As discussed above this is not what pure 
NRQCD power counting gives. We then showed that in SCET the octet structure functions 
arise naturally, and that Sudakov logarithms are summed using the renormalization group 
equations (RGEs). However, before a meaningful comparison to data can be made the color- 
singlet contribution at the endpoint must also be treated within SCET. This is the purpose 



of this work, which is an expanded version of Ref. [|1T . 

The logarithms at the endpoint for the color-singlet rate were previously studied by 
Photiadis |12| and Hautmann We find that these logarithms come entirely from coUinear 



physics, and can be summed into the form 



exp 



^a„a^(M) ln"(l 



(2) 
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Our result agrees with that of Photiadis. Hautmann argues that all logarithms cancel in 
the color-singlet contribution. We do not agree with that statement. However, the analysis 
in Ref. [|I^ was based on the eikonal approximation, which is valid for soft physics. We do 
find that there are no logarithms arising from physics at this scale, only from the collinear 
scale. 

In Section II we review the tree level spectrum, as calculated in NRQCD, and the frag- 
mentation results. In Section III we match onto SCET and present the leading operators 
consistent with the symmetries of SCET. There are two color-octet operators and one color- 
singlet operator. In Section IV we show how the color-singlet rate factorizes into hard, jet 
and usoft functions. In Section V we perform an OPE by integrating out collinear modes, 
and derive the tree-level, direct rate in the z ^ 1 limit. In Section VI we sum the Sudakov 
logarithms in the color-singlet rate, by running the color-singlet operator using the RGEs in 
SCET. By using the resummed coefficient for the color-singlet operator from the previous 
section, we have a prediction for the resummed rate in the endpoint region. In Section VII 
we discuss the phenomenology of the T radiative decay rate, including the results of the 
previous sections. We also compare our results to calculations in the literature. Finally we 
conclude in Section VIII. 



II. TREE LEVEL SPECTRUM AND FRAGMENTATION 

Before we venture to the endpoint of the photon spectrum in inclusive radiative T decay 
we review the theoretical results in the lower and intermediate range of the photon energy 
spectrum, z ^ 0.7. The leading order color-singlet contribution to the decay rate was 
first calculated in Refs. JI^. The conventional wisdom at the time was that this process 
is computable within QCD. However, Catani and Hautmann pointed out that there is 
a non-perturbative contribution which becomes important at low z. This contribution is 
due to the hadronic content of the photon, and makes itself noticed in perturbation theory 
through the presence of infrared divergences. Catani and Hautmann showed that these 
divergences can be absorbed into a non-perturbative photon structure function D{z, M). A 
consequence of their analysis is that the photon spectrum in T — X7 can be written as a 
sum of a direct contribution and a fragmentation contribution, 

dz dz dz 

where the direct term includes all contributions where the photon is produced in the hard 
scattering, and the fragmentation term is the contribution when the photon fragments from 
a parton produced in the initial hard scattering. We will look at each contribution in turn. 



A. The direct rate 

In NRQCD, the direct contribution can be calculated as an expansion in as{M), where 
M = 2mf, is the T mass, and in v, the relative velocity of the b quarks inside the bound 
state. The rate is written as 

f-^^ = J2Cn{M,z){T\On\r), (4) 
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where the C„(z, M) are short distance Wilson coefficients, calculable in perturbation theory, 
and the NRQCD matrix elements scale with a certain power in v. The lowest order contri- 
bution is the color-singlet ^Si operator. The matrix element of this operator can be related 
to the T wavefunction at the origin 

{T\O^CS,)\T) = (T|V^t^a-p£p.'T,V^p'|T) 

= §l^(0)r, (5) 
where ipl, and X-p creates a heavy quark and antiquark, respectively. 



The direct contribution was ffist calculated in the color-singlet model in Ref . , which 
for T decay is equivalent to the leading order in v in NRQCD. At lowest order in as, the 
rate is given by 

cirfA 2-z z(l-z) l-z, ,^ . 



LO 



+ 2^—\n{l - z) - 2^- ^ln(l-z), (6) 



To dz z {2-zf ^2 V ' {2-z) 



where 



32 2 ^onMMI) 

To = ^aa^e, ^ , (7) 



and Cfe = —1/3. The as correction to this rate was calculated numerically in Ref. |15|. It 
leads to small corrections over most of phase space; however in the endpoint region the 
corrections are of order the leading order (LO) contribution. 

At higher order in the velocity expansion, there are direct contributions from the color- 
octet matrix elements . The decay through a color-octet matrix element can occur at one 
lower order in as, with the bb decaying to a photon and gluon. At order v'^ and lowest order 
in as, the decay rates are 

dz ml 

— ^ — = 287refeaas j d{l-z), (9) 

dz ml 

where we have summed over J = 0, 2 in Eq. and used the relation 

{T\OsCPj)\T) = (2J+ l)(T|08(3Po)|T). (10) 

At order as, these are the only contributions. Note that these contributions are singular at 
the upper endpoint. The next-to-leading order (NLO) a^ color-octet contributions were also 
calculated in [Q. At this order, there are also contributions from the color-octet ^Si and 
^Pi channels. Large Sudakov logarithms appear in the ^5*0 and ^Po,2 contributions near the 
upper endpoint, and have been resummed in [|10|. We will discuss this resummation below. 



B. Fragmentation contribution 

Catani and Hautmann pointed out the importance of fragmentation for the photon spec- 
trum in quarkonium decays 0. The fragmentation rate can be written as 

^pfrag ^ r^dxdTa^ (11) 

dz J.. X dx '^^ \x' / ' 

o-=q,g,9 
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where the rate to produce parton a, dTa/dx, is convoluted with the probabihty that the 
parton fragments to a photon, Da-y, with energy fraction z/x. The rate to produce parton a 
can again be expanded in powers of w [^], with the leading term being the color-singlet rate 
for an T to decay to three gluons, 



LP 

dz 



^ dx dT 



X 



999 n 

dx 



M 



X 



The rate to three gluons can be obtained from Eq. 

dV 



by a change of color-factors, 



999 



dz 



5 as dVtS 
12 ela dz 



(12) 



(13) 



At order there are three color-octet fragmentation contributions, with the order par- 
tonic rates being [Q 



gg) 57raUr\Os{'So)\T) 



dz 



6{1 



dr(3pf gg) SSvra^ {T\Os('Po)\r) 



J 

dz 



mt 



6{l-z) 



dz 



mi 



6{l-z) 



(14) 
(15) 
(16) 



The as corrections to these have been calculated, and can be found in Ref. |116| . 

These rates must be convoluted with the fragmentation functions, Da-yiz, M). The M- 
dependence of the fragmentation functions can be predicted using perturbative QCD via 
Altarelli-Parisi evolution equations. However, the solution depends on non-perturbative 
fragmentation function at some input scale A, which must be measured from experiment. 



This has been done by the ALEPH collaboration for the Dq^ fragmentation function [jTT 
which is parameterized as 



Dg^{Z,fi) 



^a(jj^ 
2tt 



Pq^{z)\n 



+ C 



where C 



ln(M|/(2/iQ)) and Pq-yiz) is the Altarelli-Parisi splitting function (]T8 



Pq,{z) = 

The value for /xq extracted from the data is 



1 + (1-^) 



(17) 



(18) 



(19) 



The gluon to photon fragmentation function has not yet been measured. Here we show the 
parameterization due to Owens [|^, which uses the approximation Da-y^z, Aqcb) = 0, 



ZDg^{Z,M) 



a 

2^~ 



-0.0243(1 



ZjZ 



-0.97 



In(MVAQCD) 



2a 



0.0243('l 



«s(M)/?o 

where /3o = 11 — 2?7,//3 is the QCD beta function. 



ZjZ 



'0.97 



(20) 
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III. MATCHING ONTO SCET 



Next we turn our attention to the endpoint region. The NRQCD power-counting rules 
break down in this regime because NRQCD does not include the appropriate long distance 
modes: coUinear physics is missing from the theory. An effective theory which does include 
collinear physics is SCET ||, ^. This theory describes the interactions of highly 

energetic collinear modes with soft degrees of freedom. To describe T decay at the endpoint 
we have to couple SCET with NRQCD. 

Before we describe how to match onto SCET it is helpful to understand the scales which 
arise in the problem. Consider the momentum of a collinear particle moving near the 
lightcone. In lightcone coordinates we can write this momentum as p = {p'^ ,p~ ,p±)- Since 
the mass of the particle is much smaller than its energy, we define p"^ ~ M^A^, where M 
is the scale that sets the energy and A is a small parameter. The lightcone momentum 
components of the collinear particle are widely separated. If we choose p~ to be 0{M), then 
P_l/p~ ~ A, and p'^ /p~ ~ A^. We refer to the latter two scales as collinear and ultrasoft 
(usoft), respectively. To be concrete consider the bb pair to have momentum Mv'^ + k'^, where 
v'^ = (1, 0, 0, 0) and is (9(Aqcd) in the T center-of-mass frame. The photon momentum 
is Mzn'^/2, where we have chosen n'^ = (1, 0, 0, 1). In the endpoint region the hadronic jet 
recoiling against the photon moves in the opposite lightcone direction = (1, 0, 0, —1), with 
momentum = Mn'^ /2 + M{1 — z)n'^ /2 + . Thus the hadronic jet has n-px = Px ~ 
Next note that ^ M'^{1- z). For (1 - z) ~ t;^ ^ Aqcd/M we find 

mx ~ VMAqcd , (21) 

which is the collinear scale. This implies that for this process the collinear-soft expansion 
parameter A is of order Vl — z ~ ^JAqcd/M. The usoft scale is the component of the 
hadronic momentum in the n direction: 

2 

n ■ ~ ~ Aqcd ~ MA' . (22) 

n-px 

By matching onto SCET the large scale M is integrated out. In practice, the matching 
procedure is to calculate matrix elements in QCD, expand them in powers of A, and match 
onto products of Wilson coefficients and operators in SCET. Thus it is important to be able 
to deduce the SCET operators which can arise at a given order in A. Field theory generally 
allows any operators that are consistent with the symmetries of the theory. As explained in 
detail in Ref. 0, the symmetry of SCET which restricts the operators that can arise is the 
combined collinear- and usoft-gauge invariance of the theory. 

We will use the usoft- and collinear-gauge invariance of SCET to obtain the operators 
which are needed for the endpoint distribution of T — > X7 decays. However, we first 
review the building blocks from which these operators are constructed. In SCET there are 
fundamental fields and Wilson lines, which are built out of the fields. Furthermore there 
are two separate sectors to the theory: collinear and usoft. ^ In the collinear sector there is 
a collinear fermion field C,n,p, a collinear gluon field A!^ g, and a collinear Wilson line 



exp (-9s^n ■ An,gix)j 

perms ^ ^ 



(23) 



Soft modes do not enter our analysis so we do not include them. 



6 



The subscripts on the colhnear fields are the hghtcone direction n'^, and the large components 
of the lightcone momentum (n-g, q±). The operator projects out the momentum label [Q. 
For example n-V^n,p = 'P^n,p = n-pC,n,p- Likewise in the usoft sector there is a usoft fermion 
field Qs, a usoft gluon field A^, and a usoft Wilson line Y. Operators in SCET are constructed 
out of these objects such that they are gauge invariant. For example under coUinear-gauge 
transformations ^n,p UnC,n,p, and Wn UnWn, so the combination 

Xn ^ Wl^r^,p (24) 

is collinear-gauge invariant. This combination, however, still transforms under a usoft-gauge 
transformation Xn V"(x)xn- A collinear-gauge invariant field strength is 

^ -^W^[tV^^ + g.A^^, tV;, + 9sA^^^^,]W , (25) 
9s 

where 

= y P + n + ■ D, (26) 

and iD^^ = id^ + Qs^'^ is the usoft covariant derivative. Note that is not homogeneous 
in the power counting. The leading piece scales like A, and is given by VB^ = n^CJ^, where 
the perp subscript on B indicates that the /i index only has support over perpendicular 
components. Simplifying we obtain 

B^^ = Z!Lw\Vl + gs{A^„J^.)W. (27) 
9s 

We use these objects to build the operators we need to match onto SCET at the endpoint 
of the T — s> X'y spectrum. For further examples the reader is referred to Ref. pD[ . 

We are now in a position where we can write down the leading operators. Aside from Bj_, 
we will also need the heavy quark and antiquark fields, ipp{x) and X-pi^)^ from NRQCD. 
Consider first the color-octet ^5*0 operator |T^. The heavy quark- antiquark pair must be in 
an octet ^5*0 state, and we know that at the endpoint this pair annihilates into a photon and 
a jet. Since the operator with a quark jet has a vanishing tree level Wilson coefficient we do 
not consider it here. We only consider the gluon jet operator. So our operator must include 
the collinear-gauge invariant field strength. We want the operator to be homogeneous in 
the power counting, so it is B'^ that we use to build our operator. The most general 
LO color-octet ^5*0 operator which is consistent with these requirements, and is usoft- and 
collinear-gauge invariant is 

£prS'""H-^,/i)i?!V^P, (28) 

where Ta^i ^°\—P) is the matching coefficient. It is a function of the renormalization scale 
fi, and the large lightcone component of the colhnear momentum projected out by V. By 
momentum conservation the large lightcone momentum component must be equal to —M, 
so VB1 = -MBl, and Eq. simplifies to 

r(f;^«)(M,/x)xlp5I^p. (29) 

Note this operator is order A, and there is no operator of lower order in A which can be 
constructed. 
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+ crossed diagram 



FIG. 1: Matching the decay amplitude for bb ^ j + g in QCD and SCET. Collinear gluons are 
represented by a spring with a line through it. 

We determine the matching coefficient by requiring the SCET amphtude for a color-octet 
^5*0 bb pair decaying into a photon and colhnear gluons to be equal to the QCD amplitude 
expanded to order A. This matching is carried out order by order in an expansion in as{M). 

As an explicit example we determine Fa^^l ^°\m) to leading order in a^. The Feynman graphs 
are shown in Fig. |l]. The amplitude for a color-octet ^Sq bb pair to decay into a photon and 
a collinear gluon is simply given by the Feynman rule for the color-octet ^5*0 operator in 
Eq. (HD 

rLV'^^M) [el - ■ e)vl^T% , (30) 

where r) and ^ are two component spinors, and momentum conservation sets n - q = M. The 
tree level QCD expression for this amplitude expanded to order A is 

'^9see,ei^ (e^ - -^n ■ e) r/lpT^^p , (31) 

where = e'^'^^^fipVp. Note there is no A° piece. Setting the two expressions equal to each 
other we obtain 

F(f/^o)(M) = 2^,ee,e;^^. (32) 

The color-octet ^Pj operator is also order A. We give that operator and the matching 
coefficient in Appendix At order A these are the only two possible operators that can be 
constructed. 

Next we construct the color-singlet ^Si operator at the endpoint. A bb pair in a color- 
singlet ^5*1 configuration decays into a photon and a colhnear jet of gluons that must be 
colorless.^ The only way to construct such an operator is to include two of the B± fields in a 
colorless configuration. The only operator that can be constructed such that it is collinear- 
and usoft-gauge invariant is 

0{l, 'S,) = xlpA • a'ij^Tr{Bl T%';\V, P^) 5?} , (33) 



^ We could also have an operator where the bb pair decays to a photon and a quark- antiquark pair. The 
coefficient for this operator is zero at tree level, but would mix with the operator we have here. Numerically 
the effect of this mixing is small, and we will look into its effect in pl|. 
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FIG. 2: Matching the decay amplitude for bb ^ ^ + gg in QCD and SCET. 



where operates on the fields to the left, and A boosts from the T rest frame to a frame 
where the T has arbitrary four-momentum p2|. This operator is 0{\^) in the power count- 



ing. Once again momentum conservation forces the total momentum of the jet to be M, so 
that B1{P + P^)B^ = -MBIB^^. Introducing P_ =P-P\ Eq. (H) can be simplified to 



(34) 



The tree matching is shown in Fig. y. We obtain 



'9al3 



3M 
4:9seeb j_ 
3M ^"/^^^^ 



9,15 + -nsUf, 



(35) 



where = g^'^ — {n^n^ + n^n^^jl^ n-q^ = n-q — n-q\ and the last line holds if the photon 
is transverse. 



IV. FACTORIZATION IN INCLUSIVE T ^ X7 

Next we show that in the endpoint region the inclusive T X7 decay rate can be 
factored into a hard coefficient, a coUinear jet function, and a usoft function. We then 
perform an OPE by integrating out coUinear modes at the scale /Xc ~ Ma/1 — z ^ Aqcd, 
and match onto a usoft operator convoluted with a matching function. The usoft operator 
is nonlocal along one lightcone direction, and, as pointed out by Rothstein and Wise in 
Ref. 0, can be calculated for the color-singlet contribution. 

The factorization proof for T — > X7 is similar to that of 6 — Xg'j, and we follow 
Refs. [|], ^ . Using the optical theorem the inclusive photon energy spectrum can be written 
as 

where the forward scattering amplitude T{z) is 

T{z) = -ij d'xe-^''-^T\TJl{x)M0)\V9T ■ (37) 
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The T states are relativistically normalized, and the T indicates time ordering. We will 
show that in the endpoint region, at leading order in the SCET power counting the decay 
rate can be expressed in a factored form to all orders in as 

— = ^ H{M, li) i dk+S{k+, li) Im J^(A;+ + M(l - z), /x) . (38) 

The first step is to match the QCD current in Eq. ( |57D onto SCET operators. The 
leading order in A operators are the color-octet ^5*0 and ^Pj operators. At one order higher 
in A we match onto the color-singlet operator, Eq. (|33|) . We will show in this section how the 
color-singlet contribution factorizes, and leave the color-octet contribution for Appendix 0. 

Matching the current to leading order in A we obtain 

where the effective current is 

= ^-p^ ■ '^'^pTrji?! 5^,v- Bi], (40) 

and r[j^^^^^^(c<j) is given in Eq. (^). Note momentum conservation forces VJ{uj) = —MJ{uj), 
so we can replace V by — M in the phase factor. Substituting Eq. (|39|) into Eq. (|37D gives 

where 

and to leading order in ^^(M) we obtain 

H5S'aa'l3l3'{uJ,Uj') = - ^ 3^ '' ^ da^da'/S'Sis'- (43) 

The usoft gluons in T^^ can be decoupled from the coUinear fields as explained in Ref. |§] 
by making the field redefinition 

A'^^g = YA'^l^Y^ W = YW^^^Y\ (44) 

where the first identity implies the second. The fields with the superscript (0) do not interact 
with usoft physics. In the color-singlet contribution all usoft Wilson lines Y cancel due to 
the identity Y'^Y = 1. Furthermore the T state contains no coUinear quanta. Thus we 
can separate the coUinear physics from the usoft physics and write the effective forward 
scattering amplitude as 

x(0|TTr [Ef ^] (a:)Tr[i?f "'5^, '^'j (0)|0) , (45) 
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where we have chosen a frame such that The colhnear and usoft matrix elements 

appearing on the left hand side of the above equation can be simplified. Using rotational 
invariance 



2^, the usoft matrix element can be written as 

A' Af (T|r[^t^^^_p] (x) [£p,^^-^p.] (0)|T) = 

l^^.AfAf (T|T[^t^'=^_p] (x) [x^^^'^p'] (0)|T) . (46) 

We can then use the identity 6^^AfAj = (v^v^' — g^^'), where is the four- velocity of the T. 
Next we define a color-singlet jet function 



(OlTTr (a;)Tr[5f " ] (0)|0) = 

'^{gT'gT' + gfgr) 5^,.' j ^.^-''-'Uk^ /.) • (47) 



The jet function, J^{k^,^), is only a function of one component of the usoft momentum, 
/c"*", which follows from the colhnear Lagrangian containing only the n ■ d derivative [Q. The 
hard coefficient in Eq. (|41| ) will contract with the factors in Eqs. (|46|) and (|47|) so that the 
expression for the forward scattering amplitude can be written as 

T{z) = J2 Hiuj, fi)T,sico, z, fx) , (48) 

UJ 

where 

Teff(..,;.,/i) = y"rffc+5(x+)<5(2)(xx)et[^'^(i-^)-'="]^-^(A:+,/i) 

x(T|r[V^t^a-p](a;)[xlp,^T,V^p'](0)|T) , (49) 
and the leading order hard function is 



3 V 3M ^ 

Finally we define a color- singlet usoft function 

S{i\^^) = I ^e^^"^-(T|T[^t^.X-p](a:-)[xLp,^.^p'](0)|T), (51) 
which when substituted into Eq. (^Q]) gives the desired factored form 

Teff(o;,^,/x) = J de+J^[i+ + M{l-z),ij]S{e+,ij). (52) 
This proves the factorization theorem Eq. (|38|). 
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FIG. 3: Feynman diagram for the leading order jet function. 



V. THE OPERATOR PRODUCT EXPANSION 



As we discussed in Section III the collinear scale in T ^ X7 decay at the endpoint is of 
order ^MAqcd ~ 3 GeV. It is our opinion that this is large enough so that we can perform 
an OPE and integrate out collinear physics. However one does not have to carry out this 
step. Then only hard contributions are integrated out, and both the jet function and the 
usoft function in Eq. (P^ ) are non-perturbative. 

The OPE is carried out by expanding the jet function in powers of as{MX), and matching 
onto a non-local usoft operator of the form given in Eq. (^) convoluted with a Wilson 
coefficient. The forward scattering amplitude will then be of the form 

T{z) = j dtS{t, ^i)nj[t + M(l - /i] , (53) 

where the Wilson coefficient is labeled with a J to remind us that it contains the hard part H 
which is determined by matching SCET and QCD in an expansion in as{M). In addition it 
contains the perturbative expansion of the jet function in powers of ^^(Ma/I — z). We now 
perform the OPE for the color-singlet contribution, and finally use the results of Rothstein 
and Wise to calculate the usoft function. 

First we calculate the jet function, Eq. (|47|) , to leading order. The Feynman diagram for 
the vacuum matrix element is shown in Fig. 0. Evaluating the diagram gives 

(0|TTr[i?f-5.,P„i?fT(x)Tr[<"'5.,p_i?f^'](0)|0) = -2{gT' gf + gf 9^^' ) (54) 



(27r)4 J (27r)4 [(M + n ■ q) {K+ + k+) + ql + i5] [(M - n • q) {K+ - k+) + ql + i5] ' 

where we have used momentum conservation to set n-gi+n-g2 = and Q'fj_+Q'2_L = 0, and 
we let n-q = n-qi—n-q2 and = qi±—q2±- In addition K is the sum of residual momenta and 
k is the difference of residual momenta. In SCET there are always implicit sums over label 
momenta that come with collinear fields. This turns the integrals over residual momenta 
into integrals over the full momenta. There also is an implicit delta function which conserves 
label momenta. This delta function enforces the above momentum relations (see Ref. [§ for 
details). We get 

d'^k 1 f dk^ f dn ■ q f d'^q± 



(2^0""" 2 7 ^7 ^7 (2^- ^ ^ 



n-q q 

Once we make this replacement, the loop integral over the relative momentum in Eq. (I 
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can be performed giving 

(0|riY[sr5.,^<^(x)Tr[5r'5.,^5r'](0)|0) = '-{aT' af + gf gT' ) 



^K.^i^^ / d^JTWTJ^KT^^^e^^M^- (56) 



(27r)4 87r2 V -M^ - i6 J J_i ^ [{K+ /M){1 - ^^)Y 

By comparing the above result to Eq. (^) it is straightforward to determine the jet function, 

The decay rate is given by the imaginary part of the forward scattering amphtude. Taking 
the imaginary part of the jet function we obtain 

1 

ImJ^(A;+,/i) = — e(A;+) / d^6^,M^. (58) 

Substituting the above equation and the hard coefficient H{uj,fi) into the imaginary part of 
the forward scattering amphtude in Eq. (^51), and summing over lj gives 



2M f 1 
IxaT(z) = — I d(+ S{(+, fi) e|£+ + M{1 " ^)] / d^Hm, 

. |1 1 ,) el^ + M(l - .)] If {^^) ' . (69) 

where the 2M/M'^ accounts for the non-relativistic normahzation of the T state in the usoft 
function. The second hne is obtained by using the tree level hard coefficient Eq. (^Dp. The 
integral over ^ then gives only a numerical factor of two above. However, we will show in 
the next section that once logarithms are summed the hard coefficient depends on ^ and 
the integral is no longer trivial. Eq. (|59|) is precisely in the form given in Eq. (^31), and it is 
straightforward to read off the Hj. 

For the final step we ffist note that the usoft function may formally be written as 

Sii+,fi) = (T|^t^a-p'5(m ■ d - nxl^,cT,^p,\T) . (60) 
It is then possible to do the integral over in Eq. (pOf), giving 

ImT(z) = (T|V^t^.X-pe[m ■ g + M(l - z)]xI^,ct,^^,\T) ( ^^7?!^ ' ^^^^ 
In Ref. it was shown that 

(T|V't^,X-p0(m -d + M- 2E^)xlp,cTiijp,\T) = 

e{Mr - 2E,){T\^lcTa-pX-p'^i^p'W , (62) 

where the matrix element on the right hand side is a local operator, which can be determined 
from the T — decay rate. Remembering that z = 2E^/M we can substitute the above 
equation into Eq. (|6T|) to obtain the final expression for the imaginary part of the forward 
scattering amplitude 

167?^ 

ImT(z) = e(MT - Mz) ^Eo, (63) 

where Fq is given in Eq. (|^). Plugging the above equation into Eq. (|36|), gives the z — >• 1 
limit of the tree level rate, Eq. (Bl). 
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VI. SUMMING SUDAKOV LOGARITHMS 



One of the powers of an effective field theory approach is the abihty to sum logarithms 
using the renormalization group equations. If there exists a hierarchy of well separated scales 
in a process, then logarithms of ratios of these scales arise naturally in perturbation theory. 
If these terms are large enough so their product with the perturbative expansion parameter 
is order one, then the original expansion is no longer valid. By matching onto an effective 
theory, as we have done here, the large scale is removed, and replaced with a running scale 
/i. While the matching, which determines the Wilson coefficient, is carried out at the high 
scale, operators are run from the matching scale to the low scale using the RGEs. This sums 
all large logarithms into an overall factor, and any logarithms that arise in a perturbative 
expansion of the effective theory are order one. 

As discussed previously, near the endpoint of T — > X7 there arises a hierarchy of scales. 
As an explicit example consider the color-octet ^Sq processes. We studied the endpoint 
behavior of this contribution in a previous work [0. At next-to- leading order in ag the 
Wilson coefficient in the 2; — > 1 limit is 



1-z /, V6 3/Vl-^ 



(64) 



where C^\^So) = IQasOteln / M"^ . If this coefficient is integrated over the endpoint region 
{1 — v"^ < z < 1), the first plus distribution on the right-hand side gives rise to a double 
logarithm, log^f^, and the second plus distribution gives a single logarithm, logf^. Both 
of these are numerically of order 1/as- This clearly ruins the perturbative expansion. We 
showed in Ref. that these large logarithms are summed by matching onto SCET at the 



scale M and running to the usoft scale. The is also true of the color-octet ^Pj contribution. 

For the color-singlet ^5*1 contribution we can not explicitly give the NLO contribution 
in the endpoint region since this has only been calculated numerically [0. However, while 
the NLO correction to the color-singlet contribution does not change the LO prediction 
very much for low and intermediate values of z, it does change the endpoint prediction 
by an amount of order one. This, as we will show, is due to the presence of endpoint 
logarithms. However, unlike the color-octet case these endpoint logarithms are single, not 
double logarithms. 

In Section ^ we matched onto the SCET color-singlet operator. This integrates out the 
scale M. We now run the color-singlet operator from the hard scale to the coUinear scale, 
which sums all logarithms of 1 — 2;. Unlike the color-octet contribution the color-singlet 
operator does not run below the coUinear scale. To run the color-singlet operator given 
in Eq. (0), we calculate the counterterm for the operator, then determine the anomalous 
dimension, and finally use this in the RGEs. The four graphs that need to be evaluated to 
obtain the counterterm are shown in Fig. ^. Note that the diagrams only involve coUinear 
gluons. All of the diagrams involving usoft gluons vanish because the color-singlet process 
is transparent to usoft gluons. This is clear from the previous section where we showed that 
usoft Wilson lines completely cancel in the color-singlet contribution. This is very different 
than the color-octet contribution where usoft gluons play a crucial role. The Feynman rules 
for the vertex operators are given in Appendix 0. The result of each of the loop integrals is 
complicated so we do not give each individual expression. However, when all the terms are 
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FIG. 4: Diagrams needed to calculate the counterterm to the color-singlet operator. 



added we obtain a relatively simple result for the one-loop UV-divergent term 
e 



A 



(65) 



X 



as{fi)CA 
27r 



1 + 



M2 



M 



In h 

2M 



M 



In 



2M 



where we have introduced a sum over the index cu. Note that the divergent contribution 
depends on the large momentum component of the coUinear gluons. This means that the 
counterterm will depend on the large momentum components, as will the anomalous dimen- 
sion. This divergent piece must be canceled by Z3/Z0 — 1, where Zq is the counterterm for 
the color-singlet vertex in SCET, and Z3 is the gluon wave function counterterm 



1 + 



47r e 



This leads to 
Zo-1 



1 a. 



11 M2 + 



UJ 



M 
M + u; 



^4 - ^4 



In + 



M 



2M 



M -u 



In 



M 



UJ 



2M 



3 



(66) 



(67) 



The anomalous dimension is obtained through the standard method, and the RGE for the 
color-singlet Wilson coefficient is 



/.|^r(^'^^^)(/.,u;)=7(/^,c.)r(^'^^^)(/.,u;). 



(68) 
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We have made the dependence on the momentum label uj explicit to emphasize that the 
anomalous dimension depends on these labels. Solving the RGE gives 



In 



r(i''5i)(M,o;) 



(69) 



^^^^ 



11 M2+a;2 



M , M 
m ■ 



M + 00 



2M 



+ 



M 



M -u 



In 



M + 00 
2M 



as{M) 



where r(i''^i)(M, oo) is given in Eq. (^). Logarithms of the form ln(/i/M) have been summed 
into r*^^'^'^^)(/i, oo), and any logarithms in the operator are of the form ln(/ic//u), where fj,c ~ 
Ma/1 — z is the collinear scale. If we take ^ fic logarithms in the operator will be small, 
and all large logarithms of the ratio Hc/M will sit in the Wilson coefficient. 

We can now obtain the resummed rate, by substituting the Wilson coefficient at the 
collinear scale /i = M\/l — z, Eq. (^), into Eq. (^), giving 



lmT{z) = 2M / di+ S{i+, /i) e[£+ + M(l - z)] 



167r f4:as{M)eeb 
3M2 



(70) 



X 



dr] 



as{MVT 



as{M) 



27(r,) 



where 



2 

Jo 



11 



+ (r/2 + (l-r/)2) 



ln?7 H — ln(l — ifj 

1—7] 7] 



3 



(71) 



We have substituted ^ = 2?7 — 1. As in the previous section, using the formal result for the 
usoft function, Eq. (pOD, and the results of Ref. we obtain the final expression for the 
imaginary part of the forward scattering amplitude including resummation 



lmT{z) = Q{Mr-Mz)^Vo 



drj 



a 



■.{mVT 



as{M) 



27(r?) 



(72) 



where Tq is given in Eq. (|^). The resummed color-singlet contribution to the decay rate in 
the endpoint region is obtained by inserting the above equation into Eq. (36) 



1 dTr 



To dz 



e{Mr-Mz)z / dri 



as{M) 



(73) 



Note we can expand Eq. ([73|) in powers of as{M) and obtain an analytic expression for the 
NLO logarithmic contribution 



1 dr 

Tq dz 



^{l + [CA{2n' - 17) + 2nf] ln(l - z)| . 



(74) 



As z approaches one the 0{as) term becomes negative and of order one, precisely the 
behavior observed in Ref. |TH| for the NLO decay rate at the endpoint. 
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VII. PHENOMENOLOGY AND DISCUSSION 



In this section we combine the different contributions to obtain a prediction for the photon 
spectrum in T — > X7 decay. We will marry our expression for the color-singlet spectrum in 
the endpoint with the leading order result by interpolating between the two 

1 (iFint f 1 dV^Q \ 1 (irj-esum fnr\ 

+ ^ J— ■ (75) 



To dz \To dz J Tq dz 

The term in brackets vanishes as 2; — >■ 1, leaving only the resummed contribution in that 
region. Away from the endpoint the resummed contribution combines with the —z to give 
higher order in as{M) corrections. This clear from Eq. ([7^). 

Before we proceed we need the NRQCD matrix elements. The color-singlet matrix ele- 
ment can be related to the wavefunction at the origin, Eq. (^), which can be calculated in 
potential models, on the lattice, or extracted from the T leptonic decay rate. 



r(T ^ e+e" 



^ _ 8ajjnb)_ _ 1 Mx(n5) - 2mfe 
3 7r 3 2m/, 



'(T|Oi(35i)|T) 



ml 



(76) 



The relativistic correction was first expressed in terms of Mx — 2mb by Gremm and Kapustin 
||23|| . The matrix element is sensitive to the value of the b quark mass. We will use mj, = 4.8 
GeV, and a«(mfe) = 0.21, giving {T\Oi{^Si)\T) = 3.40 GeVl 

The color-octet matrix elements are more difficult to determine. They can, in principle, 
be extracted from data or calculated on the lattice (e.g., see Ref. |2^). NRQCD predicts that 
the color-octet matrix elements scale as compared to the color-singlet matrix element. In 
Ref. it was argued that a factor of l/2Nc should be included in a naive estimate of the 



color-octet matrix elements. For now, we leave these parameters free. 

For the color-octet ^5*0 and ^Pq rates, we need a structure function |2B[. We will use 
the simple model introduced in Ref. [l^ 



f{k+) = n(^~ ^-X^)" , (77) 

where is chosen so that the integral of the structure function is normalized to one. In 
principle the structure function can be different for the different color-octet states. But 
since we are ignorant of the non-perturbative structure function, we will naively use the 
same model for both the ^Sq and ^Pq configurations. For quarkonium the first moment of 
the structure function with respect to k'^ is 

C^l Ep,p' [V^i'^^'^r.x-p'] [x^-pT'^F.^p] |T) 

The integration limits for are from —M to Mx — M, and both A and Ai are non- 
perturbative parameters related through A = Mx — M — Ai. We use the following numbers 
in our plots: a = 1, A = 480 MeV, and Ai = -620 MeV. 

In Fig. § we show the color-singlet rate. Since each contribution is proportional to the 
same matrix element, we can normalize the rate to Fq in Eq. (|^). The dotted curve is the 
tree-level direct rate, Eq. (H), the dashed is the interpolated resummed result, Eq. ([f5|). 
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FIG. 5: The color-singlet rate. The dotted curve is the tree-level direct rate. The dashed curve is 
the interpolated resummed direct rate. The dot-dashed curve is the fragmentation contribution, and 
the solid curve is the sum of the resummed rate and the fragmentation rates. 

the dot-dashed is the fragmentation contribution, Eq. (0), and the sohd is the sum of the 
interpolated resummed and fragmentation contributions. As can be seen, the resummed 
rate rate turns over and decreases near the endpoint. At some point the resummed rate has 



a problem with the Landau pole, which can be handled in a manner described in Ref. p7 
however the Landau pole is not reached until z ~ 0.999, at which point the curve is in the 
resonance region and our results are no longer valid. 

The resummed result depends on the collinear scale. In Fig. § the collinear scale was set 
to yUc = Ma/1 — z. In Fig. ^ the solid cur ve again uses this choice of /Xc- The dotted is made 
with the collinear scale set to /Xc = Ma/2(1 — z), and the dashed has = M ^JJT^^^)J2. 
The choice of scale is a higher order effect, which can only be determined by calculating 
higher order corrections to our results. The scale variation changes the rate by about 10%. 
This can be considered part of the theoretical uncertainty. 

We also need to include the color-octet contributions to the decay rate, which requires 
knowledge of the the color-octet matrix elements. Since these are unknown at present, we 
first show in Fig. |^ the color-octet contributions with arbitrary choices of (T|(98(^>S'i)|T) 
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FIG. 6: The resummed color-singlet rate for different choices of collinear scale /ic- The solid curve 
uses He = M\/l — z, the dotted fic = M^y2{l — z), and the dashed fic = Ma/ (1 — z) /2. 



and 



{T\OsCSo)\r) + 7 



{T\Os{'Po)\r) 



mt 



(79) 



the linear combination which occurs in both the fragmentation and the resummed rates. The 
sohd curve is the octet ^Si fragmentation contribution, and the dashed curve is the octet ^5*0 
and ^Pq resummed rate convoluted with the shape function added to the octet ^Sq and ^Pq 
fragmentation rate. The ^Si octet fragmentation rate has support at larger values of z than 
the other fragmentation contributions at this order. However, higher order perturbative 
corrections to the other octet fragmentation contributions do increase the support at higher 
z 0]. The resummed octet rate produces a large peak near z ~ 0.8 which does not appear 
to be present in the data. We will therefore have to choose the combination in Eq. ( [7S| ) to 
be small enough to agree with data. 

The CLEO collaboration has measured the number of photons in inclusive T(15') radiative 



decays in Ref. p8|. The data presented does not remove the efficiency or energy resolution 
and is the number of photons in the barrel region of the detector, | cos^| < 0.7. Therefore, 
in order to compare our theoretical prediction to the data, we first integrate over the barrel 
region and then convolute with the efficiency that was modeled in the CLEO paper. We 
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FIG. 7: The color-octet rate with arbitrary choice of NRQCD matrix elements. The solid curve is 
the octet ^Si fragmentation contribution, while the dashed curve is the octet "^Sq and ^Pq resummed 
rate convoluted with the shape function and the octet ^Sq and ^Pq fragmentation rate. 

did not do a bin-to-bin smearing of our prediction. Since we do not know the size of the 
color-octet matrix elements, we will set the ^5*0 and ^Pq matrix elements to zero, and the 
^5*1 matrix element to 

{TWSr)\T) = v\T\0,{^S,)\T), (80) 

where we set v"^ = 0.08. 

In Fig. ^ we compare our prediction to the data. The error bars on the data are statistical 
only. The dashed line is the direct tree- level and fragmentation result, and the solid curve 
is the sum of the interpolated resummed result and the fragmentation result. For these two 
curves we used the value of ag extracted by CLEO from these data, as{Mr) = 0.163, which 



corresponds to as{Mz) = 0.110 ||28|. While not a perfect fit, the shape of the resummed 
result is closer to the data than the tree-level curve. We also show in this plot the interpolated 
resummed and fragmentation result, using the PDG value of as{Mz), including theoretical 
uncertainties, denoted by the shaded region. To obtain the darker band, we first varied the 
choice of nih between 4.7 GeV < ruf, < 4.9 GeV and the value of within the errors given 
in the PDG, as{Mz) = 0.1172 ±0.002 p9|. Varying nib and as modifies the extraction of the 
color-singlet matrix element in Eq. ([76|) from 3.31 GeV^ to 3.56 GeV^. We also varied the 
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FIG. 8: The inclusive radiative T photon spectrum, compared with data from CLEO [2t]. The 
error bars on the data are statistical only. The theoretical prediction, shown by the solid curve, 
includes the color-singlet interpolated resummed and fragmentation contribution and the color-octet 
^Si fragmentation contribution, with ag set to the value extracted in J^^. The other color-octet 
matrix elements are set to zero. The dashed line is the direct tree-level and fragmentation result, 
with the same choice of as- The theoretical prediction using the PDG value for ag, shown by the 
shaded bands, includes the color-singlet interpolated resummed plus fragmentation contribution and 
the color-octet ^Si fragmentation contribution. The inner band is obtained by varying m^, Og, 
and the collinear scale, while the outer band includes variation of the fragmentation function, as 
described in the paper. 

collinear scale, as we did in Fig. § from M^(l - z)/2 < fic < M^2{1 - z). Finally, the 
lighter band also includes the variation, within the errors, of the parameters for the quark 
to photon fragmentation function, Eq. (|T^), that was extracted by ALEPH [[T^. The low z 



prediction is dominated by the quark to photon fragmentation coming from the color-octet 
'^Si channel. We did not assign any error to the color-octet ^5*1 matrix element. Since it is 
unknown, there is a very large uncertainty in the lower part of the prediction that we decided 
not to show. We also do not include in our error estimate contributions from higher order 
SCET operators, which are generically suppressed by a least a power of A/M ~ \/l — z. 
Note that the color-octet ^5*0 and ^Pq contribution increases the theoretical prediction at 
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the upper endpoint, as shown in Fig. ^ It is thus clear the data favors a very small value 
for the linear combination Eq. ([79|). This is why we set these matrix elements to zero in our 
analysis. If we used small negative values for these matrix elements we would fit the shape 
in the endpoint region a bit better. 

As mentioned earlier, the NLO calculation of the photon spectrum was previously cal- 

TH]. The as correction is small for low and intermediate values 

z 



culated numerically in Ref. 



of z, but becomes larger above z ~ 0.8 until it is of order one above z ~ 0.9. Much of 
the difference between the LO and NLO piece in the endpoint is due to a aslog(l — z) 
term. However, some of the difference is due to the non-logarithmic contributions. Here 
we resummed logarithms of the form aslog(l — z) so our result includes the logarithmic 
piece of the NLO contribution. However, we do not incorporate any of the non-logarithmic 
NLO corrections. This would be done by matching onto SCET at NLO. Since some of the 
difference between the shape of the LO and NLO contribution at the endpoint is due to 
these terms, a NLO matching may give a better fit to data in the endpoint region. 

The logarithms in the color-singlet rate have previously been studied in Refs. [0, |13|. In 
Ref. |ff2l collinear logarithms were summed including mixing with quarks. Since we neglected 



mixing with quark operators (see footnote ^ we compare to the expression obtained in 
Ref. [|12| by setting the quark contribution zero. When we do that we find that the results 
are in agreement. In Ref. [|1^] it was argued that all logarithms cancel in the color-singlet 



rate, which is not what we found. However, the analysis of Ref. [|T3[ was based on the eikonal 
approximation, which is equivalent to including only usoft gluon modes. It did not include 
purely collinear modes. We do agree with the statement in Ref. |T3| that all logarithms 
due to usoft gluons cancel in the rate. This is clear in Section IV where we found that 
at leading order in A the usoft Wilson lines cancel in the usoft function. Equivalently we 
saw this cancellation in Section VI where we found that here are no diagrams renormalizing 
the color-singlet operator involving usoft gluons. However, there are logarithms due to 
collinear gluons. Note that at higher order in A, there would be diagrams with usoft gluons 
renomalizing the operator. This would correspond to logarithms in the derivative of the 
rate, whose existence was pointed out in Ref. 0]. 

In Ref. [0 a model is introduced based on the idea that the outgoing gluons shower to 
produce a jet of particles with a non-zero invariant mass. Specifically in this model gluons are 
given a mass of order rrig ~ 1 GeV, which is an estimate based on a Monte Carlo calculation 
of the invariant mass of the final state jet. It is found that this approach gives a very good 
If, as was discussed at the beginning of Section V, we treat the collinear 



fit to the data [31 



scale Ma/1 — z as non-perturbative, we would also be forced to introduce a model for the 
jet and usoft functions (since they would be non-perturbative). This approach would then 
be equivalent to the approach in Refs. |30, pT|, if we used as a model for the jet function the 



LO calculation carried out in Section V including a gluon mass of order nig ~ 1 GeV. 



VIII. CONCLUSION 

The photon energy spectrum in T X7 decay proves to be far richer than the simple 
leading order prediction based upon the color-singlet model. As was first pointed out in 
Ref. 0, at low photon energies it is crucial to include fragmentation contributions. In 
addition it has long been known that at high photon energies usoft gluon effects must be 
incorporated. In Ref. ||3^ a model was introduced to accomplish this, however, the result 



is no longer a prediction of QCD alone. In this paper we systematically treat the endpoint 
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region in the soft-collinear effective theory. Because SCET is constructed to reproduce QCD 
in the endpoint region the results of our calculation are model independent predictions of 
QCD. 

Specifically we have determined the behavior of the color-singlet contribution at the 
endpoint. This calculation consists of matching onto a color-singlet operator in SCET, and 
running the operator from the matching scale to the coUinear scale. The first step integrates 
out the hard scale set by the T mass, and the second step sums logarithms of the ratio of 
the hard and coUinear scales. Our result for the logarithms agrees with Ref. |jT^. Finally 
we obtain the decay rate by performing an OPE and matching onto a color-singlet shape 
function. As pointed out by Rothstein and Wise this shape function can be calculated. 

We combine the color-singlet contribution with the color-octet contribution calculated 
in Ref. |0 to obtain a complete theoretical prediction for the photon spectrum at the 
endpoint. Furthermore we combine our result with an NRQCD calculation that includes 
the fragmentation contribution at low photon energy to obtain a theoretical prediction for 
the entire photon spectrum. We compare our calculation to data taken by the CLEO 
collaboration p8[. Overall there is good agreement with the data. Moreover the results at 



the endpoint are now in much better agreement with data than they were previously with 
the leading order NRQCD prediction. However there is still a discrepancy between data 
and theory, and in the size and shape of the spectrum. This discrepancy might be corrected 
by including contributions that are higher order in strong coupling as and the SCET power 
counting A. 
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APPENDIX A: SCET OPERATORS AND FEYNMAN RULES 

The most general ^Pj operator is given by 

xl^T^^„s{V)BlK ■ r A ■ ^'^p , (Al) 
where the tree level Feynman rule for this operator is 

r„^.5(a;) (e^ - ■ e) A ■ p^^lpA ■ (t'T% , (A2) 

and p = p/M . To match we expand the tree level QCD amplitude to leading order in A: 

^gsee,{gi,g^^ + gi^^ + gLn'n^ - g^n^n'^) (e^ - ■ e) A ■ ^ r?lpA • ct'T% . (A3) 



n ■ q 

If we take the photon to be real {i.e. to only have _L polarization) then the n'^ term vanishes 
and the matching coefficient is 

T^^-\M) = Agsee,{g^sgi. + gLgis - ga.n'n'^) ■ (A4) 
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Next we give Feynman rules for the color-singlet operator at order g^, and order To 
obtain the g^ Feynman rule we first expand Eq. (|33D to order g"^ 

0{1, 'S,) = ^ xlpA ■ (TsiJpn ■ qn ■ q'r^'''>^{n ■q,n- q') 



{(5^ri-«,).)(^.r-«,.K)} 



n ■ An n' 



Tr<! ( -—^ql 

n ■ q ' /V n ■ q 



where Y'^^^^ is given in Eq. (pSf), which gives the following tree- level Feynman rule: 



n ■ e, 



n ■ q 

Expanding Eq. (R3|) to order g^ gives 



n ■ q 



Ogil,'S^) = igsvl^K-<Tsi^T^^'^{n-q,n-q') 

-9±-(^nJ± + 9±-g±— 



' Ann' n ■ Ann" 



n ■ q J n ■ q' n ■ {q' + q") 



which gives the Feynman rule 

-,3 



3M ^P'' 



7 7 



^ ■ (s-i + 5-2) \n- g2 n- gi 



92± 



n'n 



9i± 9L 



n- {g2 + 93) \n- gs n-g2j n ■ {gs + gi) \n ■ gi n ■ g^ 
9i± ■ 92± f 1 1 



n-gs \n ■ g2n ■ {gi + g-i) n ■ gin ■ {g2 + g^) 
92± ■ 93± 1 1 



n-gi \n ■ gsn ■ {g2 + gi) n ■ g2n ■ {g^ + gi) 
931. ■ 9i^ I 1 1 



n-g2 \n ■ gin ■ {g2 + gs) n ■ gsn ■ {g2 + gi) 



(A5) 



(A6) 



(AT) 



.(A8) 



APPENDIX B: COLOR-OCTET FACTORIZATION 

In this Appendix we show factorization for the color-octet contributions. At leading order 
in A, matching the current gives 



'(8,3Po) 



fBll 



Implicit in this formula is that we are working in a frame where the photon momentum 
defines the z axis, i.e. q± = 0. The effective theory currents are 



jatjS 



(B2) 



The Wilson coefficients are given in Eq. (|3^ ) and Eq. ([A4| ) . Next we insert Eq. (pT| ) into 
Eq. (|37D, and use momentum conservation to set the V in the phase to —M. At this order 
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in V there is no mixing between the two color-octet currents so cross-terms in the forward 
scattering amphtude vanish, and we can write 



T{z) 



(B3) 



where T'^^ is the forward scattering matrix element in the effective theory: 

T^i^so) = / rf'xe^(^^t-.)-(T|r[^;,i?2x-p'] (x) [xLp^^^p] (0) |T) g^,, (B4) 

T(t3p„) = -^| rf^xe'(*^t-^)-i(T|T[^t,i?^(p'.^)X-pO(a:) [xlp5?(p ■ (0)|T) ^i",, 

where the T states in T"^^ are non-relativistically normalized. For the sake of simplicity 
we only consider the color-octet ^5*0 contribution. The treatment of the color-octet ^Pq 
contribution is identical. The usoft gluons in T^^ can be decoupled from the coUinear fields 
using the field redefinition given in Eq. (^Ij). Under this substitution the color-octet ^5*0 
current becomes 

4^So) = Xl^ysf'-YH^ , (B5) 

where the collinear fields in Bf^ do not interact with the usoft fields. Since |T) contains no 
collinear quanta we can write the effective forward scattering amplitude as as 



^{8,iSo) 



^ j rf^xet*^(i-^)^-(T|T[^;,FT^rt^_p,](^) [^^l^FT-^yt^J (0)|T) 
x(0|TTr{T^5f "(x)}Tr{r^5f ''(0)} (7^^|0) , (B6) 
where = ^^zn^, and we used 

{T\T[i^l,YT^Yh^^,]{x) [xlprT«rVp](0)|T) = 

'-^{T\T[^l,l,YT^Y\_^]{x) [xl^YT^YH^]{m) ■ (B7) 
Next we introduce the jet function which is defined as 

(0|TTr{T^i?f "(a:)}Tr{T^i?r(0)}^;^,|0) j -^^.e-'-' J M{k% (B8) 

where the jet function Juik'^) is labeled by the large lightcone momentum of the jet (which 
in this case is M). It is a function of only one component of the usoft momentum /c"*", which 
follows from the collinear Lagrangian containing only the n ■ d derivative [Q. Inserting 
Eq. ( |B8| ) into Eq. (P^), we can integrate over three of the k momentum components to 
obtain 

x{T\T[i;l,YT^Y^X-p']ix) [xl^YT^Y^ (0)17) . (B9) 
Next we simplify the matrix element of usoft operators. We will use the following identity 

YT^Y^ = yBArpB ^ (BIO) 
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where y^"^ is the usoft Wilson hne in the adjoint representation 

r / /.X \-\AB 

37^^(0;)= Fexp iig J dsn- Af{ns)r''\ 



(Bll) 



with {T^)^^ = —if^^^. Furthermore we introduce a usoft Wilson line of finite length 

yBA^^-^yCA^Q^ = 3^^^(0, x) , (B12) 

where 



Fexp[ ig I dsn ■ A'^{ns)T^ 



EC 



(B13) 



Finally we introduce an octet usoft function defined as 



dx 



e^'^" (T| [V^;T^X-p'] y^^iO, X') [xLpT%] (0) |T) . (B14) 



Performing three of the integrals over x in Eq. ( P9|) , we can substitute the usoft function as 
defined above to obtain the desired factored form 



T-iefT 
^(8,iSo) 



di+S{i+)JM[i^ + M{1- z)] 



(B15) 



The final step is to perform an OPE by integrating out the collinear degrees of freedom. 
This is done by calculating the jet function order by order in a perturbative expansion 
in as{M^/l — z), and matching the forward scattering amplitude in SCET onto a forward 
scattering amplitude that is a convolution of a hard coefficient and a usoft operator, Eq. 
We already know the usoft operator we are matching onto, Eq. ( |B14| ). We carried out the 



matching and running for the color-octet contribution in a previous paper |T0[, and we will 
not repeat the calculation here. Readers are referred to that paper for details. 



[1] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. D 

55, 5853 (1995)] [|hep-ph/9407339i . 
[2] M. E. Luke, A. V. Manohar and I. Z. Rothstein, Phys. Rev. D 61, 074025 (2000) fep- 
I ph/9910209|] . 

[3] S. Catani and F. Hautmann, Nucl. Phys. Proc. Suppk 39BC, 359 (1995) |hep-ph/ 941 0394 1. 
[4] F. Maltoni and A. Petrelli, Phys. Rev. D 59, 074006 (1999) [|hcp-ph/9806455| . 
[5] L Z. Rothstein and M. B. Wise, Phys. Lett. B 402, 346 (1997) [|hep-ph/9701404 |. 
[6] C. W. Bauer, S. Fleming and M. Luke, Phys. Rev. D 63, 014006 (2001) [|hep-ph/0005"275| . 
[7] C. W. Bauer, S. Fleming, D. Pirjol and L W. Stewart, Phys. Rev. D 63, 114020 (2001) 
[|hep-ph/001133q| . 

[8] C. W. Bauer and I. W. Stewart, Phys. Lett. B 516, 134 (2001) |arXiv:hep-ph/0107001 |. 
[9] C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65, 054022 (2002) [|arXiv:hep- 
I ph/0109045|] . 

[10] C. W. Bauer, C. W. Chiang, S. Fleming, A. K. Leibovich and I. Low, Phys. Rev. D 64, 114014 
(2001) [|arXiv:hep-ph/0106316 |. 



26 



S. Fleming and A. K. Leibovich, arXiv:hep-ph /021 1303 . 
D. M. Photiadis, Phys. Lett. B 164, 160 (1985). 

F. Hautmann, hep-ph /010233(; . 

S. J. Brodsky, D. G. Coyne, T. A. DeGrand and R. R. Horgan, Phys. Lett. B 73, 203 (1978). 

K. Roller and T. Walsh, Nucl. Phys. B 140, 449 (1978). 

M. Kramer, Phys. Rev. D 60, 111503 (1999) |arXiv:hep-ph /9904416(| . 

A. Petrelh, M. Cacciari, M. Greco, F. Maltoni and M. L. Mangano, Nucl. Phys. B 514, 245 
(1998) [|arXiv:hep-ph/9707223|| . 

D. Buskulic et al. [ALEPH Collaboration], Z. Phys. C 69, 365 (1996). 

G. AltarelH and G. Parisi, Nucl. Phys. B 126, 298 (1977). 
J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

C. W. Bauer, S. Fleming, D. Pirjol, I. Z. Rothstein and I. W. Stewart, Phys. Rev. D 66, 

014017 (2002) iarXiv:hep-ph/0202088 |. 

S. Fleming and A. K. Leibovich, work in progress. 

E. Braaten and Y.Q. Chen, Phys. Rev. D54, 3216 (1996). 

M. Gremm and A. Kapustin, Phys. Lett. B 407, 323 (1997) [|arXiv:hep-ph/9701353 |. 

G. T. Bodwin, D. K. Sinclair and S. Kim, Phys. Rev. D 65, 054504 (2002) [|arXiv:hep 

lat/010701l[. 



A. Petrelh, M. Cacciari, M. Greco, F. Maltoni and M. L. Mangano, Nucl. Phys. B 514, 245 



(1998) [lhep-ph/9707223 |. 

S. Wolf, Phys. Rev. D 63, 074020 (2001) ||hep-ph/0010"2Tl . 

A. K. Leibovich, I. Low and I. Z. Rothstein, |hep-ph/0105066| . 

B. Nemati et al. [CLEO Collaboration], Phys. Rev. D 55, 5273 (1997) [|arXiv:hep-ex/9611020 [. 
K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 (2002). 

R. D. Field, Phys. Lett. B 133, 248 (1983). 

J. H. Field, Phys. Rev. D 66, 013013 (2002) |arXiv:hep-ph/0101158(] . 



27 



